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We introduce a new number-theoretic spin chain and explore its thermodynamics 
and connections with number theory. The energy of each spin configuration is defined in 
a translation-invariant manner in terms of the Farey fractions, and is also expressed using 
Pauli matrices. We prove that the free energy exists and exhibits a unique phase transition 

at inverse temperature (3 = 2. The free energy is the same as that of a related, non 
translation-invariant number-theoretic spin chain. Using a number-theoretic argument, 

the low-temperature ((3 > 3) state is shown to be completely magnetized for long chains. 
The number of states of energy E = log(n) summed over chain length is expressed in terms 
of a restricted divisor problem. We conjecture that its asymptotic form is (n log n), 
consistent with the phase transition at (3 = 2, and suggesting a possible connection with the 

Riemann ^-function. The spin interaction coefficients include all even many-body terms 
and are translation invariant. Computer results indicate that all the interaction 
coefficients, except the constant term, are ferromagnetic. 
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I. Introduction 

There has been considerable recent interest in the area of statistical mechanical 
models inspired by or closely connected with number theory. An overall goal of this work 
is to illuminate connections between the two disciplines, so that each might provide new 
insights and techniques useful for the other. There is, more specifically, the hope that a 
direct link between the Riemann hypothesis and the Lee-Yang theory of phase transitions, 
based on zeros of the partition function, will eventually emerge. 

In this paper, we introduce and examine a new equilibrium statistical mechanics 
spin chain model based on the number-theoretic Farey fractions. Our model is closely 
related to the number-theoretic spin chain studied by Knauf (see references given below) 
and others (see [Cv]), and in fact has the same free energy. However it differs from that 
previous work in several respects. Perhaps most important, our model is translation 
invariant by construction. In addition, our matrix formulation clarifies some of the results 
on the number theoretic spin chain. However, we have yet to conclusively demonstrate a 

connection with the Riemann ^-function, as is the case (by construction) for the previous 
model at low temperatures. 

In Section II we define the new model, and point out its connection to the number 
theoretic spin chain studied by Knauf and others, which we refer to as KSC below. We also 
exhibit exact results for the partition function at certain temperatures. Section III contains 
a proof that the free energy (of the infinite system) exists and has a unique phase transition 

at (3 = 2. In fact it is equal to the KSC free energy, establishing a direct connection between 
the two models. In Section IV we examine the number of states, and show it is related to a 
certain number-theoretic restricted divisor problem. Using a conjectured asymptotic form 
for the summed (over all chain lengths) number of states we then suggest a connection 
with the Riemann function. Section V uses the results of the previous Section to prove 
the existence of a phase transition in a different way, and also shows that the infinite 
system is in a completely magnetized state at low temperatures. In Section VI we examine 
the spin interaction coefficients, which include all even many-body terms and are 
translation invariant. Our numerical results indicate that all the interaction coefficients, 
except the constant term, are ferromagnetic. We also discuss the question of whether the 
Farey and KSC interactions are the same in the limit of long chains. 

II. Definition of the Farey Spin Chain 

We begin with a preliminary definition and then proceed to construct the Farey 
fractions. 

Definition: The mediant of the rational numbers — and — is a + c , 

J b d b+d 



Note that, for — < -and-,- e [0,11, - < iLtil < - 
b d bd 1 J b b+d d 
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Definition: The Farey fractions are defined by starting with the set F 1 



1 
l'l 



and 



recursively inserting mediants of each neighboring pair of fractions. F n denotes the 
sequence of fractions generated by this procedure up to and including the n th step; so that 



?l}' F = = {fH}' F ' = {T'H4'T Ufc ' * 



F n includes all rational fractions in [0,1] with denominators d = n, and others with 
larger denominators. Obviously F 1 cF 2 cF 3 c..cF n c F n+1 c. . We let Foo denote the 
set of all fractions generated in this manner. It follows that F„ = QJ_ [0,1]. 

Throughout this paper we employ a matrix representation of the elements of F n . 

Suppose F m = |x™,x™,...,x™ n _, }, where x™ = y, x™ m _, =^-and x™j is the immediate right 

neighbor of x™ in F m For = j = 2 m_1 , let x" 1 = — and x™. l = —, and suppose that the binary 

b (1 

expansion of j is j = (a m _,...a,a ) 2 , a k e {0,1} f or = k = m-1. Let 
a OA 



A = 



B = 



i n 

1 



= A' 



This leads to 



fc aA 



Theorem 1: M m (j) = f^A 1 " 3 ^" 1 EA 1 " Vl B a - 1 A 1 " a - 2 B V2 ...A 1 " a, B a ° = , for m * 1. 



k=0 



V d bj 



Proof. Theorem 1 follows from recursion relations for the Farey fractions very similar to 
Equation (5) of ([Ka], p. 79). 

Remark: We set a m -i = 0, so that this product always starts with A, the only exception 

( m \\ 



being when j = 2 m "! = (100...0) 2 . Then it is BA m_1 = 



employs the immediate right neighbor of 



m 



1 m-1 



m-1 1 
>1. 



In this case only, one 



1 



In what follows, we omit the i = 2 m_1 case. Thus, the fraction - is not included and 

1 

the denominator 1 only appears once in F n . Furthermore, M m (j) = AM m_1 (l), with = 1 = 
2 m " 2 ,l = (a m _ 2 ... ai a ) 2 . Therefore, for ^ j ^ 2 m "! - 1, with T m (j) defined as T m (j) = Tr(M m (j)), 
we have T m (j) = b + c = Denfx") + Num( Xj m +1 ). 
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We are now in a position to define the partition function of the spin chain. First we 
extend the definition of M m (j) to include all possible products of m factors A or B, so there 
are now 2 m matrix products. Now A = I + a., and similarly B = I + g+, where I is the (2x2) 
unit matrix and G+, a. are Pauli matrices. These satisfy o+ 2 = = a. 2 , o+o. + G-G+ = I. It 
follows that M m (j) is a linear combination of I, o+, a., o+o. and o.o+. If one exchanges A 
and B in the product, o+ and o. are exchanged, but the trace remains invariant. Thus each 
new T m (j) is exactly the same as one from the original set. 

We then interpret each T m (j) as specifying the energy of a given spin state j of a 

periodic chain of length m via E m (j) = logT m (j) . The kth spin may be regarded as down or 
up (or, equivalently, the kth site as empty or full) according to whether ak = or 1, 
respectively. 

Definition: The partition function is given by 

Z m (P) = lT m (j)- p , 

j 

where the sum extends over all 2 m matrix products, i.e. the first factor in M m (j) may be 
either A or B. 

It follows that Zi = 21-P, Z 2 = 2*"P + 2 3"P, Z 3 = 2*-P + 6 2" 2 P , Z 4 = 2*"P + 4 6"P + 8 5"P + 2 

7-P, ... 

Remark: The A-B properties of T m (j) discussed above imply that the Farey chain exhibits 
up-down (equivalently, particle-hole) symmetry. It also follows immediately from the 
properties of the trace that the energy is invariant under cyclic translation of the spin 
matrices. Therefore the spin chain has translation invariant spin interactions (see Section 

VI). 

Remark: The partition function may be evaluated exactly for certain values of (3. When (3 
= 0, Z m = 2 m , the number of states. For (3 = -1, Z m is the sum of the trace of all possible 
matrix products. Since the two operations commute, Z m (-1) = Tr( A + B) m [P]. Now A+B 
has eigenvalues 1 and 3, so Z m (-1) = l m + 3 m = 1 + 3 m . One can also calculate correlation 
functions of the A and B matrices for these (3 values in a straightforward way. These 
simplifications may also be applied to the KSC. The (3 = -1 method gives an easy way to 
derive the results of [G-K]. 

Remark: The KSC may also be expressed using the A and B matrices. To do this, one must 
replace T m (j) in the formula for the partition function with D m (j) = b = Denfx™) = (M m (j)) (2 
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and sum only over the restricted set of matrices beginning with A. We denote the 
resulting partition function Z K m -i(|3). 

In the language of [Ka], this is the canonical partition function for a chain of length 
m-1 (the leading matrix A is not counted). Since all matrix elements of M m (j) are positive, 

< D m (j) < T m (j). Thus the energy of each (restricted) state of the Farey chain is bounded 
below by the energy of the KSC. Since Z m may also be computed using the restricted sum, 
for (3 > one has Z m ((3) < 2Z K m -i((3). For (3 < 0, the inequality is reversed. 

III. Existence of the Free Energy and Phase Transition 
Definition: The Farey free energy (per spin) F(f3) is defined as 

m>= Lim ^M, 

m— >°° m 

and the KSC free energy Fk((3) is defined similarly, with Z replaced by Z K . 

Theorem 2: The Farey free energy satisfies F(f5) = Fk([3). Thus PF(f5) exists for all f5 =0,for 
(5 any negative integer, and exhibits a unique phase transition at = 2. 

Proof. For (3 = 0, (3F((3) = (3Fk((3) = -ln2 follows immediately for either spin chain by the 
remark above. For the other (3 values, we make use of Theorem 3 below, which is 

independent of the present argument. With^j- the immediate right neighbor of one 

Q 

sees that c = b + at; t = 0,1,2,3---, where the index t classifies the appearances of — with 

increasing chain length m and satisfies t < m. Here bb = l(moda) and 1 < b < a - 1 so that 
b < a < b. Hence c < (m+1) b, which we make use of as follows. First note that D m (j) = b 
and T m (j) = b + c satisfy D m (j) < T m (j) < (m + 2)D m (j). For |3 > 0, it follows that 

> Z m > ~^K^m-r Taking logarithms and dividing by -(3m gives 

(in + 2) 

fi-%,<F,„< 1 ^ + fi-i>L, 

V mj m v my 

Our results for positive (3 are then established by taking the limit m — > °°, and using the 
rigorous results that (3Fk((3) exists for all (3 > [Kd] and exhibits a unique phase transition 
at (3 = 2 [C-K]. For (3 < the same inequality holds for F m ; thus since (3Fk((3) is also known 
to exist for (3 any negative integer [C-K], the rest of the theorem is proved. 
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Remark: Another way of establishing the phase transition and gaining information on the 
behavior of the system at large |3 is given in Theorem 4. 

Remark: Note that F(-l) = ln3 by the calculation above. Making use of [C-K] gives the 
same result for Fk(-1). 

Figures 1 and 2 illustrate the free energy and energy fluctuation per spin AE (which 
is proportional to the specific heat) obtained by exact enumeration for chains up to length 
m = 16. It is clear that the convergence of the free energy with length is much slower at 

large |3, where the system is known to be completely magnetized in the infinite m limit 

(see Section V), than for small (3 values. In addition, Figure la shows that the approach to 

the limit is non-monotonic for small chain lengths, for (3=1, which is not the case at 

larger (3 values (cf. Figure lb). Finally, Figure 2 shows a peak in AE near (3 = 2, consistent 

with the phase transition there. AE increases much more rapidly with length at the peak 

((3 = 2) than at nearby (3 values. 

IV. Number of States 



In this Section we consider the summed number of states of the spin chain, that is, 
the number of states of a given energy regardless of chain length. To begin, we derive an 
expression for the number of immediate right neighbors of a given Farey fraction. 

Definition: If m is the smallest positive integer such that — e F 



then we call m the 



conductor of — . We write m = cond 



(- 
vb 







Let cond 



x x 



= m, and let — < — < — so that — ,-eF _ x are the immediate neighbors 
VbJ y' b y y' y 



of —in F m . Then we have 
b 



and 



Similarly, 



a = x + x' ; < x, x' < a 
b = y + y' ; < y, y' < b 



bx - ay = (y + y')x - (x + x')y = y'x - x'y. 



bx' - ay' = (y + y')x' - (x + x')y' = yx' - xy' 



We conclude, by induction on m, that 
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bx - ay = 1 
and bx' - ay' = -1, 

so that 



and 



bx = l(moda) ; < x < a 
ay = -l(modb) ; < y < b. 



(1) 
(2) 



Suppose b is the unique multiplicative inverse of b (mod a), i.e. the unique solution 
of(l). Then 

bb - ay = 1 so that 



y = 



bb-l 



a 



Definition: Let N t 



denote the set of all immediate right neighbors of — in F m , for m 
b J b 



cond(a/b). 

The above provides a proof of 



Theorem 3: .The immediate right neighbor of —at step m = cond 



is 



V " J 



and 



b + at . . . _ 

, t= s- t = 0,1,2, 

^bb + abt-O' 



v 



a 



Remark: If m = cond 
. a-b 



, then in F m , one can show that the immediate left neighbor of 
vb ) b 



and its immediate right neighbor is — , where a is the unique multiplicative 



is — -_ 

a b - a 

inverse of a (mod b). The form for the right neighbor is a consequence of the fact that 

aa + bb = l + ab. 



The next step is to find the number of states with a given energy, regardless of chain 
length. We proceed to count the number of solutions to 
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Tr(X) = n, (3) 

where n > 2 and X is any finite product of A's and B's beginning with A as described above. 
We count the contribution of ^ to the left side of (3) by keeping track of the 

immediate right neighbors of For this, we consider Theorem 3 and 

b + b + at = n ; t = 0,1,2,... 

which is equivalent to 

b 2 + bb + abt = bn 
or b 2 -bn + l = 0(moda). 

Therefore we look for the divisors a of bn - b 2 - 1 for 1 < a < b. There are d b (bn-b 2 -l) of 
these, where db(m) is the number of positive divisors of m that are less than b. It 
immediately follows that 

n-l 

Theorem 4: The number of solutions of (3) is <£(n) = ^d b (bn -b 2 -l). 

b=l 

Similar quantities are considered in [H-T]. 

Remark: We can put an upper bound on <£(n) by considering (4). If we ignore the 
restrictions, the number of solutions to it is just the number of unordered partitions of n 
into three positive integers. Therefore 



4>(n)< 



v 2 , 



-(n + 2)(n + l)~-n 2 , 
2 2 



where the asymptotic form applies as n — > oo. 
Conjecture: As n— >°°, 4>(n) ~ ^nlogn. 

n 

Consider the fact that ^d(m) ~ nlogn, which implies that for large n the average 

m=l 

" d „ (m) 2 _, r- 

value of d is log n, and the DDT theorem > — m sin (Vu)-n [D-D-T]. Now replace 

d(m) 71 

the quantities d(m) and db(m) by the averaged functions d(x) and d(x;y) where the latter is 
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the averaged number of divisors of y < x. Taken together, the quoted results suggest that 

f 2 _ r- 

d(x u ;x)dx sin '(Vu) ■ nlogn, and more generally the averaged quantity satisfies 



d(g;x)~-sin- 1 ( feil).logn. 

7C V log X 

log g 

for large x = n. The argument - is consistent with the fact that divisors are uniformly 

distributed on a log scale ([H-T], p. 62). The asymptotic form of the number of states then 
follows 



ii—i 

<D(n)~ Jd(b;bn-b 2 -l)db 



i 

n-1 



-MsmV 1QgP )db-logn 
71 -J \ logb(n -b) 

1 7C 

Evaluating the integral for large n, we find sin (-^=)n = — n, which leads immediately to 

a/2 4 

the conjectured formula. 

V. Thermodynamic Consequences 

The results for the (summed) number of states in Section IV may be used to derive 
some consequences for the thermodynamics. We first express the partition function as 

z rB) = 2 1 -p + 2V^4^, 



U n p 

n = 3 11 



where $ m (n) is the number of solutions of (3) for fixed chain length m, i.e. the number of 
states. The factor 2 appears in front of the summation of the preceding equation since (3) 
refers to the restricted set of matrices beginning with A. The summed number of states is 

then <D(n) = £<D m (n). If we define Z / m ((3) = Z m ((3)-2 1 " p and use the conjectured asymptotic 

m=l 

form of $(n), then 
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£z' m ((3) = 2£ 



<D(n) 



m=l n=3 11 

f, nlogn f, e(n) 2 2 " p 

n=l 11 n=l 11 J 

= -C / ((3-l) + 2e(|3)-^log2, 

where e(n) corrects the asymptotic form, e((3) is its Dirichlet transform, and £' is the 
derivative of the Riemann ^-function. If we assume that e((3) is regular, then a 
singularity in the thermodynamics for (3 = 2 follows from the pole in the Riemann 
function at (3-1 = 1. We have not been able to prove this assumption, however. 

Theorem 4: Z' m ((3 ) -> as m -> °o so that Z m (P) -> 2 2 "^, and the free energy F(f3 ) = for /? > 
3. 

Proof. We use the upper bound on <£(n) derived above. Now 

m=l n=3 11 n=3 11 

The summation in this equation is finite for (3 > 3. 

Remark: Since (3F((3) is finite at (3 = 0, Theorem 4 establishes the existence of a phase 

transition (singularity in (3F((3)) in a different way than used in Section III. In addition, 
since the partition function for long chains is given by the sum over the two lowest energy 
(j = or 2 m -l) states only, it shows that the limiting chain is in a completely magnetized 
state (all spins up or all spins down) at low temperatures. Similar behavior holds, by 

construction, for the KSC, with the partition function approaching a ratio of Riemann 
functions for (3 > 2 as m -> °°. For the Farey chain, replacing the upper bound by our 
conjectured asymptotic form for $(n), leads to F((3) = and Z m ((3) -> 2 1_ P for (3 > 2. 

VI. Spin Interaction Coefficients 

We define the spin interaction coefficients J m (t) as follows [Kb]. Let t, = t = 2 m_1 , 
specify the coefficient, and suppose that the binary expansion of t is t = (b m _ 1 ...b 1 b ) , 

m-l 

b k e{0,l}. Let j t = V^a 1 b i . Then 

i=0 

Ut) = -^rX>l) Jt E m (j), 
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where E m (j) is the energy of configuration j for a chain of length m, defined above. In this 
notation, the energy of any configuration is given via a sum over spin clusters 

E m (j) = -£(-D Jt J m (t). 

t 

Note that each factor s ; = (-l) ai e {— 1,+1} in each term may be interpreted as a spin at site i 
on the chain. The spin si is present or not in a given term according to whether bi = 1 or bi 
= 0, respectively. More explicitly, 



E m (j) = -£ 



J({b ; }). 



Each term thus defines a spin cluster, i.e. a set of sites for which bi = 1. The sites in a given 
cluster may be adjacent to one another or separated by bi = sites with no spins present. 

Lemma: The interaction coefficient J m (0) satisfies 

U0) = -^ r £lnT m (j)>-mln(^tl) = -(0. 48121.. .)m for large m. 

Proof By considering the generation of the Farey fractions F m , it is clear that each 
numerator and denominator is bounded above by the Fibonacci number F m+ i. Therefore 

/ i — \ m+1 

1 + V5 



T m (j)<2F m+1 ~-| 



, as m -> oo. 



We have calculated J m (0) by computer for chains up to m = 16 as illustrated in Figure 
3. The results indicate that it approaches -(a m + b), with a = 0.3962, for both the Farey spin 
chain and KSC. 

Since J m (0) is also (minus) the average energy for a chain of length m, we considered 
the fluctuation of the energy, i.e. 



^ 2 = ^rIanT m (j)) 2 -(^i:inT™(j)) 2 



Numerically, from calculations on chains up to length 16, a m appears to approach 0.019 m 

as m -> oo. The corresponding quantity for the KSC apparently approaches 0.014 m. Given 
the numerical uncertainties, it is possible that the asymptotic value is the same for both 
chains. 

Note that, since T m = 2, Jm(0) < 0. By the up— down symmetry discussed above, it 
follows that Jm(t) = whenever the cluster defined by t contains an odd number of spins, 

i.e. when^b; =odd, so only even interaction coefficients are non-zero. Furthermore, J m (t) 
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exhibits cyclic symmetry in t, i.e. it is invariant under translation of the spin cluster, due to 
the invariance of the trace under cyclic translation of the spin matrices mentioned above. 

Our computer results for short chains verify all the exact behavior mentioned in the 
preceding paragraph. In addition, we find several interesting features. For all non-zero 
clusters (with an even number of spins) J m (t) > 0, so all interaction coefficients (except t = 0) 
are ferromagnetic. Similar behavior occurs for the KSC in what is referred to as the grand 
canonical ensemble [Kb]. We find that as m increases, each such J m (t) apparently 
approaches a finite limit. An example is shown in Figure 4. Denoting by J(t) the limit of 

J m (t) as m -> oo, we find the values listed in Table I. Note that the (approximate) numbers 
for the pair interaction coefficients (t = (10.. .01000...) ) are consistent with a decrease of J by a 
factor of 1/2 for each increase of separation of the two spins in the cluster by one site. 
Generally, J appears to decrease with the number of spins in the cluster and the distance 

between the spins. In addition we find that in general J m (tr) ^ Jm(t)/ where t r = (b b 1 ...b m _ 1 ) 2 , 
so the interactions do not exhibit the KSC reflection symmetry. 

Theorem 2 establishes that the Farey spin chain and KSC have the same free energy. 
This raises the question as to whether J(t) is in some sense the same in both cases. This is a 
more complicated issue, in part because the KSC is not translation invariant. Our 
numerical results are consistent with the J(t) being the same, at least for KSC interactions 
with spin clusters far from the edges, i.e. with t of the form 0...0p0...0, with p = 1...1 fixed 
and each string of 0s of length proportional to m. For t of the form p0...0, so the cluster 
remains at one edge, the J(t) values are certainly different. Further, J m (t) for the KSC is 

rigorously known to exhibit the folowing behavior for m -> oo; J m (t) -> when t has an 
odd number of spins, it is small unless the length of p is small, and it has translation 
invariance [Ka]. All of these are consistent with what we find for the Farey chain, as 
described above. 

Unfortunately, the bounds relating D m (j) and T m (j) used in the proof of Theorem 2 
are not strong enough to establish equality, since each term in J m (t) will be bounded above 
and below by the corresponding term in J m K (t) plus (depending on sign) a possible term of 
magnitude ln(m+2) 2~ m , and on summation the latter gives rise to a divergent 
contribution as m -> oo . However, one can draw some conclusion about J m (0) and a m in 
this way. We assume, consistent with the Lemma and numerical results above, that both 
quantities are proportional to m in this limit. It is then easy to show that the respective 
coefficients must be the same for either spin chain, as suggested by the numerics. 

Acknowledgements. We acknowledge stimulating and useful interactions with P. 
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Table I 



t 




1 1 nnn 

1 1UUU... 


U.lJl 


1 m nnn 


n n^i ? 

U.UD1Z 


i nm nnn 


n n?Qi 


10001000... 


0.0141 


100001000... 


0.0068 


1000001000... 


0.0033 


1111000... 


0.0081 


11011000... 


0.0028 


110011000... 


0.0011 


11101000... 


0.004 
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Captions: 

Table I: Interaction coefficients. 

In 

Figure 1: Free energy F m = — vs. m for Farey spin chain at (a) (3 = 1 and (b) (3 = 4. 

m(3 

Figure 2: Fluctuation of the energy per spin 

AE^^((E m (j) 2 )-(E m (j)) 2 ) AE^-^((E m (j) 2 )-(E m (j)) 2 ) vs. (3 for the Farey spin chain with m 
= 16. 

Figure 3: J m +l(0) - Jm(0) vs. m for the Farey chain (stars) and KSC (diamonds). 
Figure 4: Pair interaction J m (11000...) vs. m for the Farey chain. 
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